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ABSTRACT 

The mean electromotive force and a effect are computed for a forced turbulent flow 
using a simple nonlinear dynamical model. The results are used to check the applica- 
bility of two basic analytic ansdtze of mean-field magnetohydrodynamics - the second 
order correlation approximation (SOCA) and the t approximation. In the numerical 
simulations the effective Reynolds number Re is 2 — 20, while the magnetic Prandtl 
number P m varies from 0.1 to 10 7 . We present evidence that the r approximation 
may be appropriate in dynamical regimes where there is a small-scale dynamo. Catas- 
trophic quenching of the a effect is found for high P m . Our results indicate that for 
high P m SOCA gives a very large value of the a coefficient compared with the "exact" 
solution. The discrepancy depends on the properties of the random force that drives 
the flow, with a larger difference occuring for 5-correlated force compared with that 
for a steady random force. 
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1 INTRODUCTION 

It is widely believed that magnetic field generation in cos- 
mic bodies is governed by turbulent motions of electrically 



conducting fluids (|Moffatt||l978| |Parker||l979| |Weiss||l994 
|Brandenburg and Su bramanian 2 005| l . One of the most im- 
portant outstanding problems of astrophysical magnetohy- 
drodynamics is to explain the phenomenon of large-scale 
magnetic activity which is observed in a wide range of astro- 
physical objects, e.g. the Sun and late-type stars, galaxies, 
accretion disks, etc. In these cases the spatial and tempo- 
ral scales of the generated magnetic fields can greatly ex- 
ceed those of the turbulent fluctuating velocity and mag- 
netic fields. According to mean-field magnetohydrodynam- 
ics ( |Moffatt|1978||Parker|1979||Krause and Radler|1980| l the 
evolution of the large-scale magnetic field B in turbulent 
highly-conducting fluid with mean velocity U is governed 
by 



dB 

~dt 



Vxf + Vx (UxB) + ?7V 2 B, 



(1) 



where the mean electromotive force, £ = (u x b) is given by 
the correlation between the fluctuating components of the 
velocity field of the plasma, u, and the fluctuating magnetic 
fields, b. We can expect a linear relationship between the 
mean electromotive force and the local large-scale magnetic 



E-mail: pip@iszf.irk.ru (VVP); mrep@cam.ac.uk (MREP) 



field, if the assumption of scale-separation holds (Moffatt 
p7^[Procto7|[2003) : 

OB 

£i = (Vx (ux b»i = a^B, +P 1] k ir L + (2) 

OXj 

where a and f3 are tensors which are usually evaluated by 
considering the dynamic equations for the small-scale veloc- 
ity and magnetic fields. If we suppose that U = 0, these 
equations are 

^ = V (u x b- (u x b) + u x B) +r;V 2 b, (3) 

i - -*-v(, + g + (f ) 

+V S Qbfc* - uu'j + i (B ■ V) b + f , 

where p is the fluctuating pressure, f is the random force 
driving the turbulence and r/, v are the molecular diffusivity 
and viscosity, respectively. 

We could also self-consistently include the effects of ro- 
tation, since the Coriolis force is linear; this enhancement is 
left for a future paper. 



It is known that the symmetric part of a and antisym- 
metric part of (3 in Q give the source and diffusion terms of 
the mean magnetic field in (JTJ, respectively. The antisym- 
metric part of a is usually interpreted as the mean pumping 
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Figure 1. The geometry of the model. 

velocity and the symmetric part of /3 may contain the ad- 
ditional source term B (e.g. Radler's, f2 x J effect, (Radler 



1969 1). For the solar dynamo the symmetric part of a (or 



simply a-effect) is a key ingredient of most mean-fields mod- 
els which claim to explain the large-scale magnetic activity 
of the Sun. 

There are currently two basic analytic methods for the 
approximate evaluation of £ and tensor coefficients in (|2j 
on the basis of (3|4 [. The most usual method is the sec- 
ond order correlation approximation (SOCA)( Krause and| 
|Radler|[l980|l which is also known as first order smoothing 
(FOSA) pIoffatt|1978] l. In this approximation, all the non- 
linear contributions of the fluctuating velocity and fluctu- 
ating magnetic fields in ( |3|4| are neglected. This approxi- 
mation has well-known limits to its accurate application. It 
is good either for poorly conducting plasma (low R m ) or 
for the weak turbulence case (low Strouhal number). Nei- 
ther limit is very appropriate in astrophysics where we have 
highly-conducting strongly turbulent fluid. On the other 
hand the r-approximation, which uses a higher order mo- 
mentum closure and could be relevant for exploring many 
common astrophysical situations, has no well defined math- 
ematically formulated limits. The particular variant of the 
r-approximation that is used in the paper will be described 
below. 

In the paper by |Courvoisier et aT| ( |2006[ > the authors 
attempted to evaluate some components of £ numerically. 
Their results indicate a nontrivial dependence of the a effect 
on the basic parameters of the turbulent flow, such as the 
correlation time, magnetic Reynolds number and the helicity 
of the flow. Here we develop a kind of shell model to explore 
some properties of mean-electromotive force and especially 
the a effect in a wide range of turbulent regimes. The model 
is useful for checking the basic approximations of mean-field 
magnetohydrodynamics as well, since it is simple enough 
to allow the rapid calculation of different cases over a wide 
parameter range while maintaining many properties of the 
full problem. 

The shell-model approach has been widely used in tur- 
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Figure 2. Electromotive force in x direction for the low(top) and 
high(bottom) P m , B/u = 0.1. 



bulence modelling ( jGledser et~aT7||1981| |Bohr et al.] [l998). 

A combination of the mean-field dynamo with a shellmodel 



was explored in Sokoloff and Prick ( 2003 l. There a dynamical 



system based on the shell-model was invoked to describe the 
dynamics of the small-scale fluctuating velocity and mag- 
netic fields. Here, we utilize a similar idea but with a dif- 
ferent purpose. Consider a velocity field with the Fourier 
representation 



u(x) = £ 



- (n) 



(k< 



+ u e 



.(k<' 



Let N = 4 and the wave- vectors form a tetrahedron as shown 
in Figure [T] Without loss of generality the wavevectors may 
be taken to have unit modulus. We suppose that the fluctu- 
ating magnetic field has the same representation, and that 
the nonlinear coupling terms only project onto this same 
set of vectors. It may be shown that the resulting closed 
nonlinear system obeys all the usual conservation laws in 
the absence of diffusion, and so seems a useful test bed for 
examining the accuracy of the various approximations. Pro- 
jecting equations (3]4l onto the given Fourier components 
we get equations for the modes: 



dth 



(i) 



— —77b 



(0 



B k 



(i) 



(5) 
(6) 



8 t u (l) = -^u (i) +l (B.k (1) )b (1) 

where the superscript means the number of the mode and 



■kj^kj. The nonlinear contributions are given in 



"ij — u *3 "i "3 

terms of the tensors A4^ and 7p" oJ\f( l > which are shown in 
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Appendix A. We suppose for simplicity that Vb = Vu = 
so that each modal eq uati on has all its terms perpendicu- 
lar to k^. Equations ( 5 6 l will be solved numerically. The 
FOSA solutions correspond to the case where all nonlinear 
contributions in (5]6l are neglected. 

To formulate the variant of the r-approximation which 
is relevant for the given model we need equations for the 
second-order products of the fluctuating fiel ds av eraged over 
the ensemble of fluctuations. Starting from (5]6l we get: 



at ( h bj 



(02,(0 



M)-'b 



+ M { ! Y b 



v- 1 (5-kW 



(0~(0 



bi u 



-(0~(0 



(7) 



(8) 
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Figure 3. Relation between the energy of the small-scale velocity 
and magnetic fields. Squares are for Case 1 and circles are for Case 
2. 
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0, ( u^bf 



(l + p-^u^f+u-lf^ 



-i f ^(0,(0 



(9) 



where the tilde above physical quantities means the complex 
conjugate and averaging over the ensemble of fluctuations is 
denoted by an overbar. In the r-approximation (see, e.g. 
Brandenburg and Subramanian (2005); Rogachevsk ii and| 
Kleeorin (2003 l ) we replace the third order contributions in 
( 7|8|9 l by the corresponding relaxation terms of the second- 
order contributions. For example, in (|9| we set 



M 



(i) AD 



+ 



r(0 



~(0 



(10) 



where r denotes the typical relaxation time of the fluctu- 
ating terms. In this formulation r is an external parame- 
ter of this approximation. We do not need to solve equa- 
tions ( 7|8|9 1. Instead we will use the left part of (10 1 to 
find the mean electromotive force obtained with the r- 
approximation. 



The model (5]6 l is clearly a good one when the diffusiv- 
ities are large, but will not give any better results than the 
other truncations when the diffusivities are small. Nonethe- 
less it does provide a useful simplification in mid ranges and 
permits the testing of the various approximations. Plainly 
a major simplification is that the fields are monochromatic. 
This could and should be remedied by increasing the number 
of shells, but this has not yet been attempted. 



THE MODEL DESIGN 



Equations (|5j |6|l were solved numerically using a second or- 
der time integration scheme. Time is measured by the typi- 





Figure 4. Correlation time. Left - hydrodynamic, right 
netic. Upper panel: Case 1; lower panel: Case 2 



mag- 



cal diffusion time, t — > t/v. The random force is normalized 
with v as well, f — > fu. 

The evolution of the small-scale velocity and magnetic 
fields depends on the typical correlation time of the random 
force. The time step is 0.003 (in dimensionless units, time is 
rescaled according to t — > ft). In what follows we consider 
two different cases. Case 1 is that of zero correlation time: 
the force is updated at each timestep. In Case 2, which has 
finite correlation time, the force was updated each 50-th time 
step. 

The effective Reynolds number is given by Re — u c i c /v. 
In computations presented below we use v — 0.05 and 
v — 0.01. We define the random driving force by writing 
f ( = w^- 1 + i\c"' x w^- 1 , and similar for initial velocity and 
magnetic fields. For each I w") is a random vector whose 
components vary between ±0.5. The term ik^- 1 x w"' is in- 
troduced to force positive helicity in the system. The initial 
velocity field is given helicity of the same sign. The electro- 
motive force associated with the (I) —mode reads 



£ 



a) 



ol (i) 



-ijq U j q 



+ C.C. 



(0 



(11) 



where tilde means the complex conjugate. Suppose the mean 
magnetic field has fixed direction, B = e x B x . The impor- 
tant component of the mean electromotive force is £ x , and 
so we define the a effect via a = £ X /B X . The mean electro- 
motive force £ is found by summation over all modes and 
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Figure 5. £ vs P m , B/v = 0.1. Left: Case 1. The FOSA solution 
is shown by circles, the exact solution by stars and squares show 
the r approximation. The right hand graph shows the same data 
for Case 2. 
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Figure 6. £ x vs B/v. The top row shows comparisons between 
the full solution of the model and the approximations for Case 1. 
At left - the case P m = 1 at right - P m = 10 3 . The FOSA solution 
is shown by circles, the full solution by stars, and squares show 
the r approximation. The bottom shows the same data for Case 
2. 



in averaging over the long-time interval equal to about 3000 
diffusion times of the system (here M is the total number of 
time-steps) 



Si 



m=M (I) 

,— T Tx (l) 



(12) 



Typical realizations of £ x = £ij q ^2 x7q f° r case Pm = 00 
and P, n = 1 are shown on the figure [2] The averaging was 
done over 16 such realizations. For the purpose of com- 
parison we also solve equations ([5] [gJ using the first or- 
der smoothing approximation (FOSA), in which the tensors 
M,N are set to zero. To test the r approximation we eval- 
uate the third order moments (see explanations above): 



(0. 
Xij 



(13) 



case v = 0.05. The results depend very much on whether 
there is a small-scale dynamo - that is whether a small scale 
field can exist in the absence of the imposed large scale 
field. The value of P m affects both the threshold and in- 
tensity of the small-scale dynamo. The relation between P m 
and the amplitude of the small-scale magnetic field fluctu- 
ations for B = is shown in Figure [3] For Case 1 there 
is a small-scale dynamo if P m ^ 10. Also, as may be veri- 
fied directly from the equations, the amplitude of the mean 
electro-motive force tends to zero if P m approaches infinity. 
This is illustrated in Figure [5]for the case B/v — 1. The typ- 
ical Reynolds number is Re m 2.2 for Case 1 and Re w 4.8 
for Case 2. 

For Case 2 the threshold is about, P m m 1. We can see 
that for large P m there is approximate equipartition between 
the energies of the fluctuating velocity and magnetic field. 

As well as examining the accuracy of FOSA, we will 
explore the usfulness of the r-approximation. The ap- 
proximation relies on knowledge of the typical relaxation 
times T (rn \T W of magnetic and hydrodynamic fluctuations. 
These quantities were found from auto-correlation functions; 
for we have 



Ah) 



(T) = 



T 



Ju) 1 ' (t)u^ (t + T)dt 
juf (t)u<'> (t)dt 
X {h) (0) 



2» (T) 



ad similarly for r 



Cm) 



r (m) 



are functions 

(h) 



First we give a detailed description of results for the 



Both t w and r 
01 B and P m , as shown in Figure [i] As we can see, r 
depends strongly on B. Its variation with P m depends on 
the existence of the small-scale dynamo. On the other hand, 
r'" 1 - 1 does not show considerable variation either with B or 
■Pm- To estimate the results of r approximation we take r = 
(tW+tW)/2. 

The dependence of the calculated mean electromotive 
force on P m for the fixed strength of B is shown in Figure^] 
The maxima are at values of P m that are close to the thresh- 
olds for the small-scale dynamo. For high P m , £ fluctuates 
strongly about zero. The dependence of the magnitude of the 
mean electromotive force on P m is not easily determined for 
small values of B because of strong fluctuations. 

To investigate the quenching of the a-effect we need 
to examine the dependence of £ x on B x . We approximate 
this with the following fitting functions, depending on three 
parameters Ai, A2, A3: 

A > B (14) 

Examples of these fits for the different cases are shown on the 
fi^6] The fit ( 14 1 does not work well for high P m as the mean 
electromotive force tends to zero and is highly fluctuating. 
However the limiting behaviour for strong magnetic fields is 
approximated quite well. The deviation of the FOSA from 
the exact solution is clearly seen for high conductivity and 
Case 2. In the same way we can say that the r approximation 
gives a better fit to the exact solution for those parameter 
values. This is confirmed by the results shown in Figure [7] 
where we show variations of ^1-3 with P m . 

Several features are quite well seen in Figure [7] First, 
in Case 1 the r approximation seems bad. Even the sign of 
effect is opposite to that for the exact solution. On the other 
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Figure 8. Plots of S x (left) and a (right) as functions of B/v 
and P m for Case 2. Case 1 is similar (sf Figur^7j. 
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Figure 9. The alpha effect (dashed line) and residual helicity 
cr( h ' (he — hie) (solid line), for the steady forcing (left) and for 
the <5-correlated random force (right), as functions of Pm find 
B/v = 0.001. 



Figure 7. Ai,2,3 vs P m - Top row: Case 1. The FOSA solution is 
shown by circles, the exact solution by stars, and the r approxi- 
mation by squares. The bottom row shows the same data for Case 
2. 



hand a significant difference between the exact solution and 
FOSA is well seen for the high P m . Second, "catastrophic 
quenching", when A2 ~ R m , is found for the high-conducting 
case. This phenomenon is more pronounced for FOSA than 
for the r approximation and the full solution. Third, in Case 
1 for FOSA the power A3 of quenching function is about 1.8 
in the whole range while for Case 2 it is slightly higher - 2. 
The quenching power of the exact and FOSA solutions are 
close. 

Plots of the amplitude of £ x and the alpha-quenching as 
functions of P m and magnetic field strength B/v are shown 
in Figure [8] Again we see that "catastrophic" quenching oc- 
curs for high P m . 

A formula that is widely quoted and has been justified 
by use of the r approximation is the simple relation between 
kinetic and current helicities in turbulent flows and the a 
effect, a ~ r (he — /ije), where he = (b • V x b) and 



h)c = (u ■ V x u)|Brandenburg and Subramanian (20051 
Krause and Radler| (1980 ); Moffatt| ( |1978| ; |Kuzanyan et al 
fl2006| l. In Figure |9| we show the a effect and residual helic- 
ity CT (h) (h c - /i;c)(with r (h) as given in Figure 4 k for two 
cases of the random force driving the turbulence. The co- 
efficient was approximately chosen to match the maximum 
magnitude of the a, we put c = 1/3 both for Case 1 and for 
Case 2. Clearly, there is no unique relation between a and 
residual helicity on the whole range of P m . Though there is 
correspondence in sign. 

Next we consider some results for a somewhat higher 



Reynolds number with v = .01. Again we present results for 
two cases. Case 1 is that of zero correlation time: the force 
is updated at each timestep and Re « 11. In Case 2, which 
has finite correlation time, the force was updated each 50-th 
time step and Re ~ 22.8. The relation between P m and the 
amplitude of the small-sca le m agnetic field fluctuations for 
B — is shown in Figure 



10 



For Case 1 there is a small- 
scale dynamo if P m 1, while it exists for P m > .1 in Case 
2 . Also in Case 2 we observe that for the high enough P m 
the energy of magnetic fluctuations is slightly larger than its 
kinetic counterpart. 

This seems to be a main reason why the a effect changes 
sign as P m varies from low to high values. Meanwhile the 
residual helicity ((he — hx,)) does not. This is demonstrated 
on Figur^TTJand helicity is shown at right side on Figur jlQ] 
The reversal of sign the a effect for high Rm was also found 
by Courvoisier et al. (2006 1. Having in mind that the en- 



ergy of magnetic fluctuations dominates the kinetic energy 
of the flow we could interpret this on the basis of results of 
analytical calculations of the a effect for a rotating strati- 
fied turbulence within r- approximation as those given in, 
e.g., Radler et al. (20031; Rogachevskii and Kleeorin (20041; 



Pipin (2007l. Suppose that the vector U characterizes the 



stratification scale and SI is a global rotation velocity then 
for the case of slowly rotating media penetrated with a weak 
large-scale magnetic field, within r- approximation we ob- 

4-7T/9 

the sign of expression in the brackets is intimately related 
to the sign of the residual helicity which is not the case for 
the computational results presented above. This point needs 
further clarifications in the multiscale model. 



tain a ~(f2 • U) 



(u 2 ^ ] . However in this theory 
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Figure 10. v = 0.01, Left: ratio between magnetic and kinetic 
energy as a function of P m - Squares are Case 1 and circles are 
Case 2. Right: the residual helicity as a function of P m - Case 1, 
dashed line, Case 2, solid line. 





Figure 11. v = 0.01, the mean-electromotive force VS Pyn with 
mean field fixed to B/v = 0.01, left is the the Case 1 right - Case 
2. We decrease the values of E obtained from FOSA to factor 10 
to make all the curves visible in one scale. Circles show FOSA, 
stars - exact solution and squares - for the r-approximation. 



DISCUSSION AND CONCLUSIONS 

One of the core issues of mean-field dynamo theory is the 
absence of reliable method for evaluation of the kinetic coef- 
ficients which describe the influence of turbulent dynamics 
on the evolution of the large-scale field. This issue is related 
to the unsolved closure problem in turbulence theories. Here 
we have attempted to construct a simple nonlinear dynami- 
cal model that can be used for this purpose. The feasibility of 
the model was demonstrated by numerical calculation of the 
nonlinear a effect. Moreover, the model is helpful to check 
two basic analytic ansatze of mean-field magnetohydrody- 
namics - SOCA(FOSA) and the r approximation, at least 
for moderate values of R m . Our results indicate that the t 
approximation may be useful in a dynamical regimes where 
the small-scale dynamo is active. On the other hand, the 
results show catastrophic quenching of the a effect for high 
P m . This is not found in analytic computations either in[Ro-| 



gachevskii and Kleeorin (20041 or in Pipin (2007 1. Certainly 



the applicability limits of this approximation need further 
clarification, but we can say with some confidence that if 



the approximation schemes fail for the present model they 
are unlikely to be very good for a fully resolved calculation. 

In the paper we present numerical calculations of the 
mean electromotive force for two different temporal regimes 
of the random force driving the turbulence. One case (Case 
1) is essentially white-noise forcing and the other (Case 2) 
is a coloured noise with a random force which was updated 
each 50th time step (for our parameters this is about two 
diffusion times of the system). We found that in the high 
conductivity limit, the difference between SOCA and the 
full solution of the model is quite significant. In particular, 
the full a effect is more than 10 times smaller than that 
from SOCA. The difference in magnetic quenching is not 

4 

very large. For high P m the aeffect is quenched a ~ B in 

— 3 

the nonlinear model though SOCA gives a ~ B which is 
consistent with previous findings by Riidiger and Kich atinov| 
fll993j );~ 



Sur et al. ( 2007 1 



The model is not competent to deal properly with tur- 
bulent diffusion because there no energy transfer to different 
spatial scales. In fact, it would be very useful to general- 
ize the simple Fourier vector space given at Figure [T] to a 
more general one with several shells. Then the effect of non- 
uniform magnetic fields and nonuniform flow on the turbu- 
lence and the mean electromotive force can be investigated 
in a similar way. One possible generalisation would be to 
consider a decomposition of the fluctuating velocity of the 
form 



»(*) = E E 



■x) 



+ CC, 



j = l,3 n=l,6 



where for the superscripts (j,n), j — 1,2,3 is related to 
the number of a vector shell and n — 1, ... 6 is the number 
of a mode, and C'C denotes complex conjugate. Each shell 
is similar to that of Figure [l] The modes of these shells 
interact in each triplet since e.g. k (M) + k (2,1) + k' 3 ' 1 ' = 
and |k (1 ' 1} | ^ |k (2a) | ^ The dynamical system thus 

obtained obeys all conservation laws. It should be suitable 
for the evaluation of a and other effects which are important 
for the mean-field dynamo, e.g., turbulent diffusion, or joint 
effect due to global rotation, nonuniform magnetic field and 
nonuniform mean flow. 
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b n bf> + bi b ( n " + b n &w + b . b t 



r(2) „ (4) _ ~(2) . (4) _ ~(3) (5) _ ~0) (5) 
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Appendix A 

Nonlinear contributions in the induction eq.(|5j 

\aW ,wfi {2) ~(3) 7 (2) ~(3) , ~( 2 )£(3) ~(2)f(3) 
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(3) 



M 



(5) 



(6) 



I (4) -(5) l (4) -(5) , ~( 4 )f(5) ~( 4 )f(5)\ 

+&„ u\> - b t u y „> +Ui b y n > -u n b\> \ , 

1.(2) (3) 1 (1) ~(3) , -(i'iO) T(l)f(3) 

i£4 ' I 6„ w> - bi u y n > +u x b y n ' - u n b\ 
+b n « t (4) - bi «i 4) + uf ) b < n > - Un ] bf ] \ , 



+ 



= ^(W'u^-bfu^+uf^-u^ 
+b n uf^ - bi + u^b^ - Unbf^ 



+ ^) fl (6)_gC3 )fi (6) +fi (3) S (6)_ fl (3) S (6) 
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